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ABSTRACT 

o . 

f^ , We study the effect that the temperature structure has on the determination of 

^N I the primordial helium abundance, Yp. We provide an equation linking r(0 iii), the 

Q ' temperature derived from the [O iii] lines, and T(He ii), the temperature of the He I 

^J ' lines, both for H II regions with O'*''^ only and for H II regions where a fraction of O"^ 

Q^ ' is present. By means of T(He ii), which is always smaller than T(0 III), we derive the 

helium abundances of 5 objects with low and very low metallicity (NGC 346, NGC 2363, 



(N 



psj . Haro 29, SBS 0335-052, and I Zw 18); these objects were selected from the literature 

^ ', because they include the 3 low metallicity objects with the best line determinations 

QQ I and the 2 objects with the lowest metallicity. Prom these abundances we obtain that 



r- 



c^ 



% 



yp(nHc) = 0.2356 ± 0.0020, a value 0.0088 lower than that derived by using r(0 iii). 



f^ ' We call this determination y^(nHc) because the collisional contribution to the Balmer 

line intensities has not been taken into account. All the recent Yp determinations in the 
_j literature have not taken into account the collisional contribution to the Balmer line 

Oj. intensities. By considering the collisional contribution to the Balmer line intensities of 

Q ! these five objects we derive that yp(+Hc)= 0.2384 ± 0.0025. 

^ ' Subject headings: galaxies: abundances — galaxies: ISM — H ii regions — ISM: abun- 

dances *** 



1. Introduction 

The determination of the pregalactic, or primordial, helium abundance by mass Yp is paramount 
for the study of cosmology, the physics of elementary particles, and the chemical evolution of galax- 
ies (e.g. Fields &: Olive 1998; Izotov et al. 1999; Peimbert &: Torres-Peimbert 1999, and references 
therein). In this paper we present a new determination of Yp based on observations of the metal- 
poor extragalactic H ii regions NGC 346, I Zw 18, NGC 2363, Haro 29, and SBS 0335-052. This 
determination is compared with those carried out by other authors. A preliminary account on the 
results for I Zw 18 and NGC 2363 is presented elsewhere (Peimbert, Peimbert, & Luridiana 2001). 

There are three problems affecting the Yp determination that need to be further analyzed: 
a) the temperature structure, b) the ionization structure, and c) the collisional excitation of the 
hydrogen lines. This paper will be mainly concerned with the effect of the temperature structure 
of the nebulae on the helium line intensities. In a future paper we will discuss in more detail 



the problem of the colhsional excitation of the hydrogen hnes (Luridiana, Peimbert, & Peimbert 
2001). The most accurate Yp values in the literature have been derived under the assumption 
of no contribution to the hydrogen Balmer lines due to colhsional excitation, we will refer to 
these determinations in this paper as y^(nHc) and to those that include the hydrogen colhsional 
effect as y^(+Hc). Due to the colhsional contribution to the Balmer line intensities the l^(nHc) 
determinations are lower limits to the 1^(+Hc) value. 

There are several pieces of evidence indicating that T(He ii) is smaller than T(0 iii): a) 
observations of NGC 346 (Peimbert, Peimbert, & Ruiz 2000, hereinafter Paper I), b) photoionization 
models of giant extragalactic H ll regions (Stasihska & Schaerer 1999; Luridiana et al. 1999, 2001; 
Relaho et al. 2001), c) in photoionization models T(0 ll) is smaller than T(0 III) for objects 
with T(0 III) higher than 12,360 K (Stasihska 1990), and a considerable fraction of the He i lines 
originates in the O"*" zone. This small difference is significative because it lowers the He^/H"^ 
determination. By combining the observed He i lines it is found that T(He ii) and iVe(He ii) are 
coupled in the sense that a lower T implies a higher iVg reducing the He^/H^ determination (see 
Paper I). 

Another problem that has to be considered to determine very accurate He/H values of a given 
H II region is its ionization structure. The total He/H value is given by: 

iV(He) _ f NeN{Re°)dV + f NeN{Re+)dV + f NeN{Re++)dV 
N{R) ~ / NeN{R'^)dV + / NeN(K+)dV ' 

rr^ z..x.„^ / N,NOie+)dV + J JV,JV(He++)dy 

= ^^^^^'^ jNjm^W ■ (^^ 

For objects of low degree of ionization it is necessary to consider the presence of He*^ inside 
the H"*" zone, while for objects of high degree of ionization it is necessary to consider the possible 
presence of H" inside the He"*" zone. For objects of low degree of ionization /Ci^(He) might be 
larger than 1.00, while for objects of high degree of ionization /C-F(He) might be smaller than 
1.00. This problem has been discussed by many authors (e.g. Shields 1974; Stasihska 1983; Peha 
1986; Vilchez & Pagel 1988; Pagel et al. 1992; Armour et al. 1999; Peimbert & Peimbert 2000; 
Viegas, Gruenwald, & Steigman 2000; Viegas & Gruenwald 2000; Ballantyne, Ferland, & Martin 
2000; Sauer & Jedamzik 2001). 

Each H II region is different and a good photoionization model is needed to estimate an 
/C-F(He) of very high accuracy. In this paper we will assume that the H II regions are chemically 
homogeneous. 

In section 2 we will adapt the t^ and Tq formalism introduced by Peimbert (1967) to relate 
T(0 III) and T(0 ii) with r(He ii), the full mathematical treatment is presented in the Appendix. 
In section 3 we will rediscuss the Y determination of NGC 346, the most luminous H II region in 
the SMC, carried out in Paper I. In sections 4 to 7 we redetermine the He"''/H+ values for I Zw 18, 



NGC 2363, Haro 29, and SBS 0335-052, based on the observations of Izotov, Thuan, & Lipovetsky 
(1997) and Izotov et al. (1999) and on photoionization models for these objects. These four objects 
were selected from their sample based on the following criteria: the two objects with the lowest 
heavy element abundances and the two objects with the smallest observational errors (Peimbert & 
Peimbert 2001). We derive He+/H+, A''e(He ii), r(He ii), and t(3889) self-consistently based on 
all the observed He i line intensities, with the exception of those strongly affected by underlying 
absorption. By combining the He^/H"*" values of these four objects with that of NGC 346 we derive 
a new l^(nHc) in section 8. The collisional contribution to the Balmer lines is estimated in section 
9. The discussion and conclusions are presented in sections 10 and 11. 



2. Temperatures 

2.1. T(0 III) and r(He ii) for a pure 0++ nebula 

We first consider an H II region where the He ii, the O iii, and the H ii volumes coincide 
perfectly. If we further assume that the temperature is constant we have a very simple model 
which gives us a first approximation to a real nebula. 

If we now assume that the H II region has an average temperature given by 



/ TNeNp dV 



and a mean square temperature variation given by: 

2 ^ J{T-nYN,N^dV 

T^jN.NpdV ' ^^' 

we find that the observed r(4363/5007), derived from the 7(4363) //(5007) ratio, is given by (Pe- 
imbert 1967): 

^.^ N ^ f [90800 
r(0 III) = To 1 + 



T^ ^ -)> (4) 

and that the temperature associated to the He i recombination lines, that originate in the He ii 
region, is (see Paper I): 

r(He II) = To (1 - 1.4^2) . (5) 

This means that the temperature that should be used to measure the helium abundance is given 
by: 

nHe„,^.(0.„,(:-[^ 



0.2 



(6) 



In Figure 1, and based on equation 6, we present T(He ii)/T(0 iii) versus t^ for different values 
of T(0 III). 



2.2. T(He ii) for a nebula with 0+ and 0++ 

We consider now another H ii region where the He ii volume coincides with the H ii volume, 
but now oxygen is present both in the form of O^^ and O^, with the O^ weight described by the 
fraction a: 

^ jN,N{0+)dV 

J NeN{0+) dV + f NeN{0++) dV ^ ' 



If the average temperatures in the O ii and O iii zones are given by: 

jTNeN{0+)dV 



!02 



j NeN{0+)dV 

and 



(8) 



_ JTN,N{0++)dV 
°^ fNeN{0++)dV ^' 

respectively, then the average temperature for the whole ionized region Tq is equal to 

To = aTo2 + (l-a)To3; (10) 

the fractional difference between Tqs and Tq can thus be characterized by the parameter [3 so that 

^ = 1 + A (11) 

and the ratio between Tq2 and Tq is: 

^ = 1-^(1^)=! + ,,-^. (12) 

In the general case when the O ii and O iii zones are not of uniform temperature, but include 
temperature fluctuations, represented by: 

2_ /(T-To2)^JVeJV(0+)dy 

2 T^^jNeN{0+)dV ^ ' 



and 



. _ /(T - ro3)2iVeiV(0++) dV 



*3 T^J-N,N{0++)dV ^^^^ 



respectively, the mean square temperature fluctuation for the entire H ii region, t , is given by: 

,2 _ M02 - Jo \ ,^ nMoS — J0\ , ,2 M02 \ , /, n .2 M03 



were the first two terms represent the deviation from uniform temperature given by the difference 
in temperature between each zone and the average for the whole H ii region, while the other two 
represent the additional contribution due to inhomogeneities in each zone. 



Substituting the first 2 terms with equations 12 and 11 and by keeping the last 2 terms to 
second order we obtain: 

2 



af3 



P' 



a 
1 — Q 



■^ ' " ^ + (1 - a) /?2 + Qt2 + (1 - a) tl 



a 



+ ati + (l-a)t 
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Tq3 and t| are related to T(0 iii) by (see equation 4): 



r(o III) = To3 1 + 



90800 



Tf 



03 



(16) 



(17) 



while Tq2 and t^ are related to an O ii temperature, determined from the /(3727)//(7325) ratio, 
that in the low density limit is given by (Peimbert 1967): 



T(0 II) 



f02 



1 + 



97300 



Tf 



02 



(18) 



It is convenient to obtain an equation, with the same physical meaning than that of equations 
17 and 18, that relates Tq and t^ with a temperature that can be derived from observations, 
particularly with those of the O II and O III forbidden lines. Intuitively we can define the following 
equation (the formal derivation of equation 19 is given in the Appendix): 

t2- 



T(0 ii+iii) = To 1 + 



n 



t' 



(19) 



where Tq and i^ correspond to the values for the entire volume, 6 represents an effective average 
over the excitation energy needed to produce the oxygen lines, 



e = (1- a) 90800 + a 97300 



(20) 



and T(0 ii+iii) is an average between T(0 ii) and T(0 iii). However, it can be shown (see Ap- 
pendix) that, the simple average of the temperatures weighted by a and 1 — a does not give an 
equation accurate to second order and that T(0 ll+lll) needs a factor F of the form: 



r(0 ii+iii) = [(1 - a) T(0 iii) + a T{0 ii)] T, 



(21) 



r is very close to unity for objects where q(1 — q) < 0.08 or where T(0 ii) is similar to T(0 ii) 
(0.9 < r(0 ii)/T(0 III) < 1.1); but for H ii regions where there are appreciable fractions of both 
O"'" and O^'^ and where there is an appreciable difference between T(0 III) and T(0 ii), F must 
be calculated as described in the Appendix. (All computations done in this paper were done using 
the formalism described in the Appendix.) 



Finally, solving equation 19 for Tq and substituting it in equation 5 it follows that 

^2^ 



r(He ii) = T(0 ii+iii) 1 



T(0 ii+iii) 



0.2 



i" 



(22) 



which can now be used to determine the He abundances. To derive r(0 ii+iii) from the observed 
T(0 ii) and T(0 iii) values we need to use equations 7, 20, and 21 for T = 1; for T larger than one 
we also need to use equations A2, A3, and A7. If only T(0 iii) is available we need a photoionization 
model to obtain a relationship between T(He ii) and T(0 iii). 

In Figure 2 we present the relationship between T(He ii) and T(0 iii) as a function of t^ for 
objects with different O"*" fractions. Figure 2 is based on equations 7, 20, 21, A2, A3, and A7 and 
the relationship 

To2 = 2430 + ro3(1.031 - To3/54350), (23) 

derived from the photoionization models by Stasiiiska (1990). From this equation it is obtained that: 
for To3 = 12,360 K then Tq^ = Tq2, for higher temperatures Tqs > T02, and for lower temperatures 

To3 < Tq2. 

From Figure 2 and equation 23 we note the following cases for objects with T(0 iii) > 12,360 K: 
a) For Q = 0.00 and ts^ = 0.000, then r(He 11) = r(0 iii) and t^ = o.OOO. b) For a / 0.00 and 
t2^ = ^3^ = 0.000, then T(IIe 11) is always smaller than T(0 iii) and t^ is always larger than 0.000. 
c) For a = 0.00, and ta^ > 0.000, then r(He 11) < re(0 iii) and t^ > o.OOO. d) For a / 0.00 and 
^2^ ^ ^3^ > 0.000, then T(IIe ii) is even smaller than in case b and t^ is even larger than in case b. 

Very metal-poor extragalactic H 11 regions, those that have been used to derive the primordial 
helium abundance, have T(0 III) larger than 12,360 K and are best represented by case d. 

Figure 2 is also related to the Averaged Quantities block of Cloudy's output (Ferland 1996; 
Ferland et al. 1998; Kingdon & Ferland 1995). In that section of Cloudy's output there is a line 
labeled "Peimbert" with the following entries: "T(0 Illr)" that corresponds to T(0 iii) in this 
paper; "T(Bac)" is the hydrogen temperature resulting from the predicted Balmer jump and H/? 
for the optically thick case, corresponds to the outer surface of the object and does not represent 
typical H 11 regions; "T(IIth)" is the hydrogen temperature resulting from the predicted Balmer 
jump and H/? for the optically thin case, corresponds to T(Bac) defined by Peimbert (1967) and used 
by many authors, and applies to typical H 11 regions; "t^(Hstr)" is the structural value computed 
across the H 11 region, and it is called t^ in this paper. "T(03-BAC)" and "t2(03-BC)" are the 
average temperature and the mean square temperature fluctuation values derived from the predicted 
T(0 III) and T(Bac) values, they correspond to Tq and t^(Hstr) only when all the O is in the 0++ 
stage. In the presence of O"*" and according to Figure 2 and equation 23 we can distinguish two 
cases: when r(0 iii) is higher than 12,360 K then "t2(03-BC)" is higher than t^, but when r(0 iii) 
is smaller than 12,360 K then "t2(03-BC)" is smaller than t^. 

From photoionization models of H 11 regions it has been found that t^ is in the 0.002 to 0.03 
range, with typical values around 0.005 (e.g. Gruenwald & Viegas 1992; Kingdon & Ferland 1995; 
Perez 1997); from observations it has been found that t^ is in the 0.01 to 0.04 range, with typical 
values around 0.03 (e.g. Peimbert 1995, 1999, and references therein). As a first approximation 
due to equation 23 we expect minimum t^ values from photoionized models with T ~ 12,360 K 
and an increase in t^ for hotter and cooler models. Most of the observational determinations of t^ 



come from H II regions with T(0 III) ~ 8000 K which are biased toward one of the higher ends 
of the t^ distribution. Very metal-poor H ii regions are biased to the other high end of the t^ 
distribution, particularly those with WR stars and a large fraction of O^. The t^ values derived 
from photoionization models correspond to T(He ii) values 3% to 11% smaller than r(0 iii). In 
the presence of additional sources of energy, to those provided by photoionization, the difference 
between r(He ii) and T(0 iii) might be larger. 



3. NGC 346 

In Paper I we have analyzed NGC 346, the most luminous H ii region in the Small Magellanic 
Cloud. Due to the relatively small distance to NGC 346 it is possible to place the observing slit 
avoiding the brightest stars, thus minimizing the effect of the underlying absorption in the IIe''~/II"'" 
determination. 

The determination of Yp based on NGC 346 can have at least four significant advantages and 
one disadvantage with respect to those based on distant H ii region complexes: a) no underlying 
absorption correction for the helium lines is needed because the ionizing stars can be excluded from 
the observing slit, h) the determination of the helium ionization correction factor can be estimated 
by observing different lines of sight, c) the accuracy of the determination can be estimated by 
comparing the results derived from different lines of sight, d) the electron temperature is generally 
smaller than those of metal poorer H ii regions reducing the effect of coUisional excitation from the 
metastable 2^ S level of He I and from the ground level of H i, and e) the disadvantage is that the 
correction due to the chemical evolution of the SMC is in general larger than for the other systems. 

From the observations of the He I lines AA 3889, 4026, 4387, 4471, 4922, 5876, 6678, 7065, 
and 7281 (see Table 1) and from the temperatures associated with the O III lines and the Balmer 
discontinuity in region A we have used a maximum likelihood method to determine simultaneously 
and self-consistently several parameters. Our results are: A'^g(He ii)= 144]^3g cm~^, T(He ii)= 
11,860 ± 370 K, He/H = 0.07910 ± 0.00072, and r(3889) = 0.0. For this temperature range we 
obtain that 0/H= (141 ib 21) x 10^^. In Table 2 we present He/H and x^ as a function of several 
A''e(He ii) and T(He ii) values, which include a few representative temperatures and the densities 
favored by those temperatures. From this table it can be seen that the minimum x^ values imply 
a strong correlation between the temperature and the density, in the sense that the lower the 
temperature, the higher the density, and consequently the lower the He/H value. In Paper I it is 
also found that Te(He ll) is 8.6% it 3% smaller than Te(0 III). The parameters that represent this 
model are presented in Table 3. 

In Table 4 we present the y(nHc) values for different temperatures and densities, where we 
have included the He^^/H^ ratio (see Paper I) and an /CF(He) = 1.00. 

Relaho et al. (2001) have computed photoionization models for NGC 346 with Cloudy, they 
have found: a) that there is a negligible amount of H*^ inside the He"*" zone, in agreement with 



Paper I, b) that T(He ii) is 5% smaller than T(0 iii) in their best model, and c) that the T(0 iii) 
computed by the models is 9% smaller than the observed value, indicating the presence of additional 
heating sources not included in the photoionization models and probably implying that the real 
difference between T(He ll) and T(0 lll) is higher than predicted by the models. 



4. I Zw 18 

Up to now, with the exception of NGC 346, it has not been possible to derive the /C-F(He), 
A^e(He ii), T(He ii), r(3889), and the He"'"/H+ ratio based only on the helium lines. Usually the 
observed T(0 iii) value has been used to complement the information provided by the He i lines, 
moreover it has also been assumed that the observed T(0 iii) is equal to T(He ii). For I Zw 18, 
SBS 0335-052, NGC 2363, and Haro 29 we will make use of the observed T(0 iii), but we will not 
assume that it is equal to T(He ii), instead we will make use of photoionization models computed 
with Cloudy to relate both quantities. 

In general photoionization models predict T(0 iii) values smaller than observed (Stasiiiska & 
Schaerer 1999; Luridiana et al. 1999; Luridiana & Peimbert 2001; Luridiana et al. 2001; Relaho 
et al. 2001) indicating the possible presence of an additional heating source not considered by the 
models; this result also implies that the T(He ii) values predicted by the models will be smaller than 
those T(He ii) values derived from the observed T(0 iii) values. Since the lower the temperature, 
the lower the derived He/H value (see for example Table 5), the use of the temperatures predicted 
by the models yields spuriously low He/H values. On the other hand the models can give us an 
estimate of t^(0 iii), probably a lower limit. The formalism presented in section 2, together with 
i^(0 III), provides a relationship between the observed T(0 III) and r(He ii). 

From photoionization models of I Zw 18 based on Cloudy (Luridiana et al. 2001) we find that: 
/CF(He) = 1.00, r(3889) = 0.010 ±0.005, and ^^(O iii)= 0.013 ±0.006. We consider t^{0 iii) as a 
lower limit to t^ (see equations 16 and A9). From the formalism of section 2 together with a = 0.265 
and T(0 iii) = 19, 060 K (Izotov et al. 1999) we obtain t^ = 0.024 ± 0.006 and consequently that 
r(He ii) is 10.2% smaller than T(0 iii). 

To derive the He'*"/H+ ratio we have used the maximum likelihood method, MLM (see Paper I). 
The inputs are: a) r(3889) = 0.010 ± 0.005, b) t^ = 0.024 ± 0.006, c) r(0 iii)= 19, 060 K and d) 
the observations of AA 3889, 4026, 4471, 5876, 6678, and 7065 by Izotov et al. (1999), where we 
have adopted the line intensity ratios presented in Table 1, corrected by underlying absorptions of 
3.3A, 0.8A, and 0.8A for AA 3889, 4026 and 4471 respectively, that correspond to a burst model of 
4 Myr (Gonzalez-Delgado, Leitherer, & Heckman 1999). The underlying absorptions predicted by 
the burst model have been reduced by about 10% to take into account the presence of the nebular 
continuum emission; a similar reduction has been applied to the underlying absorptions predicted 
by the starburst models for NGC 2363, Haro 29, and SBS 0335-052. 

To obtain He+/H+ values we need a set of effective recombination coefficients for the He and H 



lines, the contribution due to collisional excitation to the hehum line intensities, and an estimate of 
the optical depth effects for the helium lines. The recombination coefficients that we used were those 
by Storey & Hummer (1995) for H, and Smits (1996) for He. The He i collisional contribution was 
estimated from Kingdon & Ferland (1995) and Benjamin, Skillman, & Smits (1999). The optical 
depth effects in the He i triplet lines were estimated from the computations by Robbins (1968), 
and based on Paper I we will assume that the He i singlet lines are produced under case B. 

The maximum likelihood solution amounts to a) Ne(iie ll)= ST^^g cm~'^, b) He"'"/H"'" = 
0.07673 ±0.00312, c) x^ = 4.2, and d) Te(Re ii)= 17, 120 ±600 K, 10.2% smaller than r(0 iii) (see 
Table 3). For this temperature range we obtain that 0/H = (19 ± 2) x 10^^. There is not enough 
information to independently derive the t^ and r(3889) values from the observed lines in this object, 
therefore their output values are equal to the input values. Table 5 presents the y(nHc) values for 
four t^ values: 0.024 (the optimum value), 0.030 (± la), 0.018 (- la), and 0.006 (the minimum 
possible value given by the size and temperature of the O^ zone); and 6 different densities that 
include those provided by the MLM for each t^. The y(nHc) values include the contribution due 
to He++ and an ICF(Re) = 1.00. 



5. NGC 2363 

The T(0 III) values for the four objects observed by Izotov et al. (1997) and Izotov et al. 
(1999) were derived by us from the observed line intensities, we obtain the same results for I Zw 18 
and SBS 0335-052 than Izotov et al. (1999), but we obtain values 650 K higher for NGC 2363 and 
Haro 29 than Izotov et al. (1997). 

Luridiana et al. (1999) produced detailed photoionization models of NGC 2363, they find also 
that the T(0 ill) predicted by the models is considerably smaller than 15,750 K, the observed value; 
from their models they find also that t^ is in the 0.007 to 0.029 range. This amounts to a T(He ii) 
6% ± 2% smaller than r(0 ill). From the models of NGC 2363 and the slit used by Izotov et al. 
(1997) we also find that: ICF(He) = 0.993, indicating the presence of neutral hydrogen inside the 
ionized helium region (see equation 1), and r(3889) = 0.50 ± 0.15. 

Support for the relative high density comes from the Ne[S ii] derived by Izotov et al. (1997) 
and by Esteban et al. (2001) that amount to 120 cm^^ and 340 ± 120 cm~^ respectively and the N^, 
[Ar iv] values in the 400 to 800 cm~^ range derived by Perez, Gonzalez-Delgado, & Vilchez (2001) 
and the 1300 ± 500 cm~^ value derived by Esteban et al. (2001); the [S ii] lines trace the outer 
parts while the [Ar iv] lines trace the inner parts of the H ll region. 

For the MLM we use the following inputs: a) r(3889) = 0.5 ± 0.3, b) t^ = 0.018 ± 0.011, c) 
r(0 ill)= 15,750 K, and d) the observations of AA 3820, 3889, 4026, 4387, 4471, 5876, 6678, 7065, 
and 7281 by Izotov et al. (1997) slightly modified to take into account underlying absorptions of 
0.3A, 1.8A, 0.4A, O.2A, and 0.4A for AA 3820, 3889, 4026, 4387, and 4471 respectively (see Table 1), 
that correspond to a burst model of 3 Myr (Gonzalez-Delgado et al. 1999). We did not consider 



- 10- 

A4922 because its intensity was many a away from the solution increasing x^ to unacceptable 
values; we also adopted an error for A3889 three times larger than the one given by Izotov et al. 
(1997) due to our estimate of the combined errors produced by the reddening correction and the 
underlying absorption correction. 

With these inputs we obtain a) Ne{Yie ii)= 250+^^ cm^^, h) t^ = 0.021 ± 0.010, c) t(3889) = 
0.98 ± 0.25, d) He+/H+ = 0.07997 ± 0.00153, e) x^ = 12-9, and /) T(He ii) = 14, 710 ± 440 K, a 
value 6.5% smaller than T(0 iii) (see Table 3). 

Note that the t^ and r(3889) output values are not the same as the input values. This is 
because the quality of the measurements for NGC 2363 is better than that of I Zw 18, and thus 
the values of t^ and r(3889) can be directly estimated from the lines, however these values are not 
as accurate as those of NGC 346 and the error bars would be very large, therefore the observed 
line intensities get mixed with the t^ and t(3889) input values in the MLM, giving us a different 
output. 

For this temperature distribution we obtain 0/H= (95 ± 8) x 10~^. Table 6 presents the 
y(nHc) values for four t^ values: 0.021 (the optimum value), 0.031 (+ la), 0.011 (- Icr), and 0.001 
(the minimum possible value given by the size and temperature of the O'*' zone); and 6 different 
densities that include those provided by the MLM for each t^. The y(nHc) values include the 
contribution due to He"^'^ and an /CF(He) = 0.993. 



6. Haro 29 

Haro 29, also known as I ZW 36 or Markarian 209, is one of the two blue compact galaxies 
with the smallest observational errors in the list of Izotov et al. (1997), the other being NGC 2363. 

From photoionization models of Haro 29 based on Cloudy by Luridiana et al. (2001) and 
the slit used by Izotov et al. (1997) it is again found that the T(0 iii) values predicted by the 
models are in the 14,550 to 15,100 K range, values considerably smaller than the observed value 
of 16,050 K; from these models it is also found that ICF(He) = 0.995 - 0.996, r(3889) ^ 1.4, and 
t^ = 0.002 -0.004. 

For the MLM we use the following inputs: a) r(3889) = 1.43 ±0.715, h) t^ = 0.020 ±0.007, and 
c) the observations of AA 3820, 3889, 4026, 4387, 4471, 5876, 6678, 7065, and 7281 by Izotov et al. 
(1997) slightly modified to take into account underlying absorptions of 0.4A, 2.2A, 0.6A, 0.3A, and 
0.6A for AA 3820, 3889, 4026, 4387, and 4471 respectively (see Table 1), that correspond to a burst 
model of 3.2 Myr (Gonzalez-Delgado et al. 1999). We did not consider A 4922 because its intensity 
was many a away from the solution increasing x^ to unacceptable values; we also adopted errors 
for AA3889 and 7065 three and two times larger than the one given by Izotov et al. (1997) due 
to our estimate of the combined errors produced by the reddening correction and the underlying 
absorption correction. 
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From these inputs we obtain: iVe(He ii)= 45lgg cm"^^ r(3889) = 0.899 ± 0.296, t^ = 0.019 ± 
0.007, and He"'"/H+ = 0.08157 it 0.00188. The density result is of great concern, not only is the 
lower limit unphysical, but the derived value is similar to the root mean square density, Ne{rms), 
that amounts to = 47 it 7 cm~^. Giant H ii regions usually show N(,{rms) values considerably 
smaller than those densities determined through a forbidden-line ratio or from the He I lines based 
on the maximum likelihood method, N^ {local); to quantify the difference the filling factor e is used 
and is defined by means of the relation: 

N^{rms) = eN^ (local). (24) 

Typical values of e for giant H ii regions are in the 0.001 to 0.1 range (e.g. Paper I; Luridiana et al. 
1999; Luridiana & Peimbert 2001; Relaho et al. 2001) that would imply a A''e(He ii) in the 149 to 
1490 cm~^ range. Since the density derived from the MLM implies an e ~ 1.00, an unacceptable 
value, we decided to make another estimate of A''e(He ii). 

From the Cloudy models by Luridiana et al. (2001) we found that a good fit to the observed 
line intensities is obtained with e ~ 0.04, which corresponds to an A''e(He ii)= 235 it 85 cm~^. 
Fixing this density as our preferred value we obtain r(3889) = 0.41 it 0.30, t^ = 0.023 it 0.007 (see 
Table 3), and He^/H+ = 0.07800 it 0.00178, other values can be obtained or interpolated from 
Table 7. The y(nHc) values include the contribution due to He"'""'' and an ICF(He) = 0.995. 



-6 



For this temperature distribution we obtain 0/H = (78 it 10) x 10 



7. SBS 0335 052 

After I Zw 18, SBS 0335-052 is the extragalactic H ii region with the second lowest metallicity 
known; and being brighter than I Zw 18 makes it critical in determining Yp. 

We have computed a series of Cloudy models for SBS 0335-052 (Luridiana et al. 2001). We 
find that the T(0 III) predicted by the models is T(0 iii) = 17,500 - 18,500 K, again considerably 
smaller than the observed value of 20,500 K. For the three central values of the observing slit used 
by Izotov et al. (1999) (center, 0."6SW, 0."6NE) we have extracted from the models the following 
predictions: ICF(He) = 0.9985, r(3889) ^ 1.26, and t^ in the 0.010 to 0.015 range. 

To correct for underlying absorption the blue helium lines and the Balmer lines observed by 
Izotov et al. (1999) we adopted the values predicted by Gonzalez-Delgado et al. (1999) for a 3 Myr 
old instantaneous burst with Z = 0.001. The corrections for AA 4026, 4471, 4922, 3889 and H/3 
amount to 0.4A, 0.4A 0.3A 1.8A and 2.0A respectively. Note that Izotov et al. (1999) corrected 
the Balmer lines for an underlying absorption of 0.2 A and that they did not correct the helium 
lines. We adopted the C(H/3) value determined by Izotov et al. (1999) that is mainly based on the 
/(Ha)//(H/3) ratio, this ratio is not affected by underlying absorption (but it might be affected by 
coUisional excitation of the Balmer lines, see below). The line ratios are presented in Table 1. 
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From the MLM and the fohowing inputs: a r(3889) = 1.26 ±0.63, b) t^ = 0.020 ±0.007, and c) 
the modified observations of AA 3889, 4026, 4471, 5876, 6678, and 7065 (A 4922 was not considered 
because its intensity was many a away from the solution increasing the x^ to unacceptable values); 
we obtain iVe(He ii)= 297+^^ cm'^, He+/H+ = 0.07640 ± 0.00152, r(3889) = 1.60 ± 0.35, and 
f = 0.021 ± 0.007, this t^ corresponds to a T(He ii) = 19, 290 ± 310 K (see Table 3). 

For this temperature distribution we obtain 0/H= (24 ± 3) x 10^®. Table 8 presents the 
y(nHc) values for four t^ values: 0.021 (the optimum value), 0.028 (± Icj), 0.014 (- Icr), and 0.004 
(the minimum possible value given by the size and temperature of the O"*" zone); and 6 different 
densities that include those provided by the MLM for each t^. The y(nHc) values include the 
contribution due to He"*"^ and an IC-F(He) = 0.9985. 



8. Determination of y^(nHc) 

In Table 9 we compare the y(nHc) maximum likelihood values computed in this paper for 
t^ / 0.000 (see Tables 4-8) with the y(nHc) values for t^ = 0.000 computed by Izotov et al. (1997, 
1999). The differences in the Y values between both types of determinations are of the order of 
0.0070, and are manly due to the treatment of the temperature. The errors that we quote are larger 
than those by Izotov and colaborators, both groups include the errors in the observed line intensity 
ratios, however we also include the uncertanties on the A''e(He ii), t^, and r(3889) values while 
they do not. None of the y(nllc) errors include the uncertanties due to the collisional excitation 
of the Balmer lines (which will be discussed in the next section) nor the prescence of other possible 
systematic effects. In Figure 3 we also present the y(nHc) versus 0/H diagram for all the objects. 

To determine the l^(nllc) value for all the objects it is necessary to estimate the fraction of 
helium present in the interstellar medium produced by galactic chemical evolution. We will assume 
that 

^^ = ^"^A0' ^^^^ 

where O is the oxygen abundance by mass. The AO baseline given by the five objects in the sample 
is very small and consequently produces large errors in the Ay/AO determination, therefore we 
decided to adopt Ay/AO = 3.5 ± 0.9, the slope derived in Paper I. From this slope and the five 
y(nHc) values for t^ / 0.000 presented in Table 9 we derive yp(nHc)= 0.2356 ± 0.0020. In Table 9 
we present the l^(nHc) determination derived from the y(nHc) values for t^ = 0.000 obtained by 
Izotov & Thuan (1998) and Izotov et al. (1999) for Haro 29, NGC 2363, I Zw 18, and SBS 0335- 
052, and find that our result is 0.0088 smaller than theirs. The main difference is due to our use of 
r(He ll) instead of r(0 iii) used by them. 
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9. Collisional excitation of the hydrogen Unes 

Davidson & Kinman (1985) were the first to estimate the colHsional contribution to the Balmer 
hnes and its effect on the determination of Yp] they made a crude estimate of this effect for 
I Zw 18 and concluded that the contribution to /(Ha) may be roughly 2%. All the subsequent 
determinations of Yp in the literature have been derived under the assumption of no contribution 
to the hydrogen Balmer lines due to collisional excitation, we have referred to these determinations 
in this paper as y^(nHc) and as y^(+Hc) to those primordial helium abundance determinations 
that include the collisional contribution effect. Further discussion of the relevance of the collisional 
contribution to the Balmer lines is presented elsewhere (Stasihska & Izotov 2001; Luridiana et al. 
2001). 

From observations it is not easy to estimate this effect for two reasons: small changes are 
expected in the Balmer line ratios, and the increase in /(Ha)//(H/3) due to collisions can be ascribed 
to a spuriously higher C(H/3). The He/H line ratios are affected by the increase of the Balmer line 
intensities due to collisions and by the decrease of C(H/3) due to the increase of /(Ha)//(H/3). 

From a series of Cloudy models it is found that the collisional contribution to /(H/3) for 
I Zw 18 and SBS 0335-052 is in the 2% to 6% range, for Haro 29 and NGC 2363 in the 1% to 
3% range and for NGC 346 in the 0.6% to 1.2% range. The effect of collisions on the Balmer 
line intensities together with the Cloudy models for I Zw 18, SBS 0335-052, and Haro 29 will be 
discussed extensively elsewhere (Luridiana et al. 2001). The Cloudy models for NGC 2363 and 
NGC 346 are those by Luridiana et al. (1999) and Relaiio et al. (2001) respectively. 

To estimate y(+Hc) for NGC 346, NGC 2363, Haro 29, and I Zw 18 we used the published 
Balmer line intensities without modifying the C(H/3) determination, the reason is that the published 
C(H/9) values are already very low and do not seem to indicate that they have been overestimated. 
For NGC 346 the published C(H/3) value is 0.15ib0.01 and the value derived from the embedded 
stellar cluster amounts to 0.19 it 0.01 (see Paper I), for NGC 2363 and Haro 29 C(H/3) amounts to 
0.11 and 0.00 respectively (Izotov et al. 1997), and for I Zw 18 it amounts to 0.015 ± 0.020 (Izotov 
et al. 1999). The y(+Hc) estimates are presented in Table 9, where the collisional contribution 
was obtained from the models mentioned in the previous paragraph. The errors in the y(+Hc) 
determinations include the same errors as the y(nHc) determinations plus the error produced by 
the uncertainty in the collisional contribution effect, but they do not include the errors due to other 
possible systematic effects. 

Alternatively Izotov et al. (1999) derived for the three central positions of SBS 0335-052 a 
C(H/3) value of 0.237 it 0.019, but the value derived from /(Ha)//(H/3) is higher than those derived 
from I (Hj) / 1 (Hf3) and /(H5)//(H/3) (after correcting the Balmer lines for underlying absorption 
as mentioned in section 7); this difference could be due to collisional excitation of the Balmer 
lines, therefore from the corrected Balmer lines for underlying absorption and from the effect of 
collisions, as predicted by one of our Cloudy models, we obtained a C(H/3) of 0.15; for this value 
of the reddening we corrected all the line intensities and from them we derived the y(+Hc) value 
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presented in Table 9. 

Our preliminary result of 1^(+Hc), is about 0.0028 larger than l^(nHc) (see Table 9). From 
Table 9 it can be seen that the y(+Hc) values for NGC 346, NGC 2363 and Haro 29 are similar 
to those for I Zw 18 and SBS 0335-052. This result may indicate that the corrections adopted for 
coUisional excitation of the Balmer lines have been overestimated because the production of helium 
due to galactic chemical evolution is expected to be higher for the first three objects than for the 
last two (see equation 25). 



10. Discussion 

The brightest extragalactic H ii regions with temperatures in the 14,000 to 16,000 K range 
might be the best objects to determine Yp. The reasons are the following: a) there are many of 
these objects available and the brightest of them are brighter than the metal-poorest objects known 
which implies that the errors in the line intensity ratios will be smaller, and b) the effect of collisional 
excitation of the Balmer lines is considerably smaller than for the metal-poorer objects (those with 
temperatures in the 18,000 to 22,000 K range) reducing the error in the Y determination due to 
this effect. 

On the other hand the correction to the Y determination to obtain Yp based on objects with 
temperatures in the 14,000 to 16,000 K range increases the error in the Yp determination, but 
in general the size of this error is expected to be smaller than those introduced in the previous 
paragraph because the correction is small and can be obtained with high accuracy under the 
assumption that AY/AO = 3.5 it 0.9 (see Paper I). From chemical evolution models it is found 
that changes in AY/AO are small even in the presence of a very strong burst of star formation 
(Carigi, Colin, & Peimbert 1999); moreover these objects seem to have been forming stars over 
periods of at least 1 to 2 Gyr implying that most of their oxygen formed before the present burst, 
as in the case of Haro 29 (e.g. Schulte-Ladbeck et al. 2001), and consequently that the expected 
changes in AY/AO due to bursts are very small. 

The y^(nHc) value derived by us is significantly smaller than the value derived by Izotov & 
Thuan (1998), from the Y - 0/H linear regression for a sample of 45 BCGs, and by Izotov et al. 
(1999), from the average for the two most metal deficient galaxies known (I Zw 18 and SBS 0335- 
052), that amount to 0.2443 ± 0.0015 and 0.2452 ± 0.0015 respectively. 

Our Te(He ll) determinations, that we prefer, produce Y values that are about 0.007 smaller 
than those derived by Izotov & Thuan (1998) and Izotov et al. (1999), a small but significative 
difference, as we will see below. It should be noted that the abundances for NGC 2363, SBS 0335- 
052, Haro 29, and I Zw 18 obtained by both groups are based on the same observations. 

From photoionization models computed with Cloudy we estimate that T(He ii), the temper- 
ature that should be used to determine the helium abundance, should be at least 5% smaller than 
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T(0 III). Moreover, if in addition to photoionization there is additional energy injected to the H ii 
region the difference between T(0 iii) and T(He ii) might be even larger. 

Figure 4 shows the helium, deuterium, and lithium abundances predicted by standard Big 
Bang nucleosynthesis computations with three light neutrino species for different values of rj, the 
baryon to photon ratio (Thomas et al. 1994; Fiorentini et al. 1998), also in this figure we present 
observational abundance determinations of these elements. The implications of this figure for the 
determination of $7^, the baryon content of the Universe, are presented in Table 10. 



11. Conclusions 

Based on the best observations of extragalactic H ii regions of low metallicity available in the 
literature we redetermine the primordial helium abundance, l^(nHc), taking into account all the 
observed He i line intensities with the exception of those strongly affected by underlying absorption. 
We derive He/H, iVe(He ii), r(He ii), and t(3889) self-consistently. 

From the same data for NGC 2363, SBS 0335-052, Haro 29, and I Zw 18 we obtain y(nHc) 
values about 0.007 smaller than those derived by Izotov, Thuan and collaborators, the differences 
are small but significant and are mainly due to the lower T(He ii) values adopted by us. In the 
self-consistent solutions the lower T(He ii) values imply higher densities; the higher the density, the 
higher the collisional contribution to the He I line intensities and consequently the lower the helium 
abundances (see Tables 2 and 4-8). The y^(nHc) value derived by us from NGC 346, NGC 2363, 
SBS 0335-052, Haro 29, and I Zw 18 is in good agreement with that derived from NGC 346 in 
Paper I. 

The y(nHc) values are lower limits to the real Y values because they do not consider the 
excitation of the Balmer lines due to collisions. A preliminary estimate of y(+Hc) is presented in 
Table 9. In a future paper we will give a full discussion of this problem (Luridiana et al. 2001) 
based on detailed models for each object. 

The y^(nHc) of 0.2356 it 0.0020(1cj) combined with standard Big Bang nucleosynthesis com- 
putations (Thomas et al. 1994; Fiorentini et al. 1998) implies that, at the 2a confidence level, Of, 
is in the 0.011 to 0.020 range for h = 0.7 (see Table 10). Also in Table 10 we present seven other 
determinations of ilf, not based on Yp, our result is in good agreement with five of them, those 
based on: the primordial Lithium determination. Lip, the Cosmic Background Imager, the baryon 
budget in the galactic vicinity (z = 0), the baryon budget at z ~ 3.00 and that derived under 
the assumption of a flat universe, the SNIa magnitude redshift data to derive ^m, and a baryon 
fraction estimate based on X-ray observations of rich clusters of galaxies (this determination is 
labeled SNIa in Table 10); on the other hand our result is in disagreement with two of them, those 
based on: the primordial deuterium determination, Dp, and BOOMERANG. 

Our 1^(+Hc) value of 0.2384 it 0.0025 is also in good agreement with the five determinations 
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of rift mentioned above and is closer than the y^(nHc) value to the other two determinations but 
still in disagreement with them. 

The errors in all determinations need to be reduced to constrain even further the domain 
available to non standard BBN models. 

It is a pleasure to acknowledge several fruitful discussions with: G. Ferland, B. E. J. Pagel, E. 
Skillman, G. Stasihska, G. Steigman, S. Torres- Peimbert, and S. Viegas. 



A. Second order approximation to the two temperature scheme 

In this appendix we prove the validity of equation 19. To do this, we need to give a physical 
meaning to T(0 ll+lll) and consequently to T, presented in equation 21. Therefore we need to 
derive the functional form of T that will make equation 19 accurate to second order. 

Equation 4 is designed (and by analogy 17 and 18) to take into account the increase in emissivity 
of the forbidden O iii lines with increasing temperature, therefore giving higher T(0 iii) — Tq 
differences to regions with larger temperature fluctuations. We want to give the same physical 
meaning to equation 19, thus we want the T(0 ii+iii) — Tq difference to increase with t^. 

This can be physically explained by examining equation 16 where the three terms are positive 
and thus increase t^. The first term is related to the temperature difference between the O ii and 
O III zones and the last two terms are related to the temperature fluctuations in each zone. The 
fluctuations associated with the last two terms correspond to the r(0 ii) — T02 and T(0 iii) — Tqs 
Gw3o3h)W differences, therefore by adopting a r(0 ii+iii) of the form (a)T(O 11) + (1 - a)T(0 iii) 
the fluctuations will produce a T(0 ii+iii) — Tq difference. The first term should also contribute 
to the T(0 ii+iii) — Tq difference, but since this term does not represent a change in t^ or t|, 
T(0 11) and T(0 lll) will not change either, consequently another term is needed in the definition 
of T(0 ii+iii) This effect is introduced by multiplying the direct average by a factor T, which 
represents the second order correction due to the difference in temperature between the O 11 and 
O III zones, thus T(0 ll+lll) will be of the form: 

r(0 ii+iii) = [(1 - a) r(0 III) + aTiO 11)] T, (Al) 

where 

r = 1 + d^-f- (A2) 

Here S is the relative difference between T(0 11) and T(0 iii) and 7 represents a (free) parameter 
which will permit to recover the physical meaning of equations 17 and 18. Thus S is given by 

« = ?^-l- (A3) 

r(o 11) ^ ' 
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To have a second order correction in /3 it is enough to determine 5 to first order; to first order 
r(0 III) = Tos and r(0 ii) = T02, therefore 



P 



a 



(A4) 



Substituting equations 17 and 18 into equation Al we obtain 



T(0 ii+iii) 



:i-a)(l + /3)(l-3|j+(l-a)^| 



+ a{ l + P 



r 



P 



a 



1_3( (l_a)^ + a^ 



f2 

2 



1-3^) +a 

2 



97300 q 



+ — (1 - q) 90800^ + a97300^ 
To V 2 2 



(A5) 



substituting T and keeping this equation accurate to second order in the temperature variations 
(/9) *2, H: and t) we obtain 



r(0 ii+iii) 
To 



1 + 5^ 
-3((l-a)ti + at2)l 

+7^ ((1 - a) 90800t| + a97300ti) -. 
Tq 2 



(A6) 



To simphfy this equation and recover the physical meaning of equation 4 it is necessary for 7 



to be 



7 = -(l — a)a-- — 3, 
' 2 ^ V^o ' 



(A7) 



here it is enough to derive 7 to zeroth order, so dividing 9 by any one of the available temperatures 
(To, T02, To3, r(0 11) or T(0 iii)) would be approximately right, thus we obtain 



T(0 ii+iii) /, ,n a ^^\-a\l 



To 



+- 



90800 



a J 2 



a 2 



97300 / 2 ^2l-aA 1 



(A8) 



Under the assumption that 



/2 
^2 



''3> 



(A9) 



equation 16 becomes 
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t^ = tl + (i^L^^ (AlO) 

a 



and finally equation A8 becomes 

T(0 ii+iii) = To ( 1 + ;^-3 -), (All) 



1 + 


[tA 


t2^ 
2y 


). 


e 


-3 


t2 


T(0 ii+m) 


2 



or 

To = r(0 ii+iii) ( 1 - ——^ - 3 - ) . (A12) 

^ ^ V L^(0 ii+iii) J 2 y ^ ^ 

Equation A12 can now be combined with those equations that relate Tq and t^, to some observable 
temperature like r(Bac) or T(He ii). 

To adopt a value of F = 1 in equation 21 is equivalent to assume that the first term in 
equation 16 is equal to zero. We have already mentioned the parameter space where it is a good 
approximation to assume F = 1. To obtain very accurate He/H values the exact value of F, that 
includes the second order effects implied by t^, should be used. 
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Table 1. Adopted He i line intensities relative to H/3 



He I line 


NGC 346 


I Zw 18 


NGC 2363 


Haro 29 


SBS 0335 - 052 


3820 






0.012±0.00r 


o.oioio.oor 




3889^ 


0.0940±0.0017^ 


0.0898±0.0071^ 


0.087±0.00r 


0.095±0.00r 


0.0924±0.0028^'-'' 


4026 


0.0185±0.0007^ 


0.0202±0.0036^ 


0.018±0.00r 


0.020±0.00r 


0.0173±0.0006^'^ 


4387 


0.0047±0.0002^ 




0.005±0.00r 


0.005±0.001^ 




4471 


0.0384±0.0005^ 


0.0387±0.0025^ 


0.040±0.001= 


0.039±0.001^ 


0.0382±0.0007'=-^ 


4922 


0.0100±0.0002^ 




0.013±0.00r 


e 


0.0093±0.0004^'^ 


5876 


0.1064±0.0012'' 


0.0936±0.0028<^ 


0.106±0.001'^ 


0.102±0.001<^ 


0.1027±0.0017^ 


6678 


0.0296±0.0002'' 


0.0263±0.0018<^ 


0.029±0.001'^ 


0.029±0.001<^ 


0.0263±0.0006^ 


7065 


0.0211±0.0002'^ 


0.0239±0.0016^ 


0.029±0.001'i 


0.024±0.001'i 


0.0358±0.0007f 


7281 


0.0063±0.0003^ 




0.006±0.001'i 


0.005±0.001'i 





'^Where the contribution due to the H8 line has been subtracted. 

^From Peimbert et al. (2000). 

'^From Izotov et al. (1997, 1999), but modified due to underlying absorption, see text. 

■^From Izotov et al. (1997, 1999). 

'^Presented by Izotov et al. (1997), but not used, see text. 

Modified by assuming a different H/3 underlying absorption from that used by Izotov et al. 
(1999), see text. 
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Table 2. 7V(He+)/iV(H+)^ and x^ for Region A 

iVe(He ll) 
(cm~^) 
r(He II) 

t^ (K) 69 90 110 144 182 220 

0.001^ 12,800 8052 8014 7980 7924 7865 7810 

(12.4)^ (14.0) (18.4) (31.7) (54.2) (83.3) 

0.014 12,240 8036 8001 7969 7917 7861 7809 

(13.0) (8.81) (7.53)^= (10.8) (21.7) (39.0) 

0.022 11,860 8026 7992 7961 7911 7858 7808 

(22.3) (14.5) (9.58) (6.53)'=''^ (10.0) (19.8) 

0.030 11,480 8015 7983 7953 7906 7855 7807 

(39.0) (27.9) (19.7) (10.8) (7.53)= (10.3) 

^Given in units of 10^^, for the case of no hydrogen colhsions, x^ values 
in parenthesis. 

'^Minimum t^ value, assuming that ^2^ = ^3^ = 0.000 and considering the 
difference between T02 and T03, see text. 

=The minimum x? value at a given t^. 

The smallest x^ value for all t^'s and densities, thus defining T(He 11) 
and A^e(He 11). 
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Table 3. Parameters of the preferred model for each H II region. 



Parameter 


NGC 346 


IZw 18 


NGC 2363 


Haro 29 


SBS 0335 - 052 


r(0 iii)^ 


13070 ± 50 


19060 ± 610 


15750 ± 100 


16050 ± 100 


20500 ± 200 


r(He ii)^ 


11860 ± 370 


17120 ± 600 


14710 ± 440 


14880 ± 300 


19290 ± 310 


^2b 


0.022 ± 0.008 


0.024 ± 0.006 


0.021 ±0.010 


0.023 ± 0.007 


0.021 ±0.007 


7V(He ii)^ 


144t^^ 


87+94 


250t7j 


235 ± 85^ 


29711 


r(3889)^ 


0.000 


0.010 ±0.005 


0.98 ±0.25 


0.41 ± 0.30 


1.60 ±0.35 


ICF{Rey 


1.000 


1.000 


0.993 


0.995 


0.9985 



'^Observed value. 

•^Output from the MLM fit, see text. 

^Prom photoionization models, see text. 
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Table 4. y(nHc) values for NGC 346 











iVe(He II) 














(cm-3) 








T(He il) 

(K) 












t2 


69 


90 


110 144 


182 


220 


0.001^ 


12,800 


0.2436^ 


0.2428 


0.2420 0.2407 


0.2393 


0.2381 


0.014 


12,240 


0.2433 


0.2425 


0.2417^ 0.2405 


0.2392 


0.2380 


0.022 


11,860 


0.2430 


0.2423 


0.2415 0.2404^' 


■=0.2392 


0.2380 


0.030 


11,480 


0.2428 


0.2420 


0.2414 0.2403 


0.2391^ 


0.2380 



'^Minimum t^ value, assuming that t2^ = ^3^ = 0.000 and considering the 
difference between Tq2 and T03, see text. 

''Minimum x^ value at a given t^, see Table 2. 
■^The smallest x^ value for this object, see Table 2. 
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Table 5. y(nHc) values for I Zw 18 



iVe(He II) 

,-3\ 



cm 



T(He II) 



t^ 


(K) 


0.006^ 


17,920 


0.018 


17,390 


0.024 


17,120 


0.030 


16,860 



25 59 77 87 97 150 

0.2404 0.2375'" 0.2362 0.2354 0.2346 0.2307 

0.2400 0.2373 0.2359*' 0.2351 0.2344 0.2307 

0.2398 0.2371 0.2358 0.2350^'^= 0.2343 0.2307 

0.2396 0.2370 0.2357 0.2349 0.2342'^ 0.2307 



'^Minimum t^ value, assuming that ^2^ = ^3^ = 0.000 and considering the 
difference between Tq2 and T03, see text. 

^Minimum x^ value at a given t^. 

■^Best fit (minimum x^ value) for this object. 
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Table 6. y(nHc) values for NGC 2363 



iVe(He II) 

,-3\ 



cm 



T(He II) 

t^ (K) 100 183 214 250 289 338 

O.OOP 15,620 0.2500 0.2459'" 0.2445 0.2430 0.2414 0.2395 

0.011 15,160 0.2493 0.2455 0.2442'" 0.2427 0.2412 0.2394 

0.021 14,710 0.2487 0.2451 0.2438 0.2424^''=0.2410 0.2393 

0.031 14,260 0.2480 0.2446 0.2434 0.2421 0.2407'" 0.2391 

'^Minimum t^ value, assuming that ^2^ = ^3^ = 0.000 and considering the 
difference between Tq2 and T03, see text. 

^Minimum x^ value at a given t^. 

■^Best fit (minimum x^ value) for this object. 
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Table 7. y(nHc) values for Haro 29 











Ne{B 


:e ll) 














(cm 


-') 








r(He II) 

(K) 














t2 


50 


100 


150 


235 


320 


500 


O.OOP 


15,860 


0.2477 


0.2451 


0.2427 


0.2390 


0.2356 


0.2296 


0.016 


15,190 


0.2471 


0.2446 


0.2424 


0.2389 


0.2356 


0.2200 


0.023 


14,880 


0.2468 


0.2444 


0.2422 


0.23881^ 


0.2356 


0.2302 


0.030 


14,570 


0.2465 


0.2441 


0.2421 


0.2388 


0.2356 


0.2304 



'^Minimum t^ value, assuming that t2^ = ^3^ = 0.000 and considering the 
difference between Tq2 and T03, see text. 

^Preferred value, see text. 
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Table 8. y(nHc) values for SBS 0335-052 



iVe(He II) 
(cm^^) 



r(He II) 



t2 


(K) 


0.004^ 


20,020 


0.014 


19,590 


0.021 


19,290 


0.028 


18,990 



241 257 279 297 316 362 

0.2439 0.2425^ 0.2409 0.2396 0.2382 0.2351 

0.2442 0.2430 0.2413^ 0.2400 0.2387 0.2356 

0.2445 0.2433 0.2416 0.2403^''= 0.2390 0.2360 

0.2447 0.2435 0.2419 0.2406 0.2393^ 0.2363 



'^Minimum t^ value, assuming that ^2^ = ^3^ = 0.000 and considering the 
difference between Tq2 and T03, see text. 

^Minimum x^ value at a given t^. 

■^Best fit (minimum x^ value) for this object. 
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Table 9. Y Comparison 







y(nHc)^ 


y(nHc)^ 


y(+Hc)= 




Qa 


f = 0.000 


t^ / 0.000 


t^ / 0.000 


NGC 346 


1714 ± 255 




0.2405 ± 0.0017 


0.2421 ± 0.0018 


NGC 2363 


1155 ± 97 


0.2456 ± 0.0008 


0.2424 ± 0.0035 


0.2456 ± 0.0039 


Haro 29 


948 ± 122 


0.2509 ± 0.0012 


0.2388 ± 0.0040 


0.2420 ± 0.0044 


SBS 0335-052 


392 ± 36 


0.2463 ± 0.0015 


0.2403 ± 0.0048 


0.2450 ± 0.0060 


IZwl8 


243 ± 24 


0.2429 ± 0.0070 


0.2350 ± 0.0072 


0.2412 ± 0.0085 


1^ (sample) 




0.2445 ± 0.0009^^ 


0.2356 ± 0.0020<^ 


0.2386 ± 0.0025"^ 



'^Oxygen abundance by mass for t^ rsj 0.02, in units of 10 ^. 
^Izotov & Thuan (1998); Izotov et al. (1999). 
'^This paper. 

Derived from the sample under the assumption that AY/AO = 3.5 ± 0.9 (see 
Paper I). 
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Table 10. Baryon densities for h = 0.70 



Method 




n 


b 




Source 


Dp 




0.035 - 


■ 0.049 


{2a) 


O'Meara et al. (2001) 


Lip 




0.013 - 


■ 0.026 


{2a) 


Suzuki et al. (2000) 


Boomerang 




0.031 - 


■ 0.059 


{2a) 


Netterfield et al. (2001) 


CBI 




0.006 - 


• 0.060 


(niin-niax) 


Padin et al. (2001) 


(z = 0) 




0.007 - 


■ 0.041 


(min-max) 


Fukugita et al. (1998) 


{z ~ 3.0) 




0.010 - 


■ 0.060 


(min-max) 


Fukugita et al. (1998) 


SNIa 




0.018 - 


■ 0.066 


{2a) 


Steigman et al. (2001) 


Fp(nHc),t2 = 


0.000 


0.025 - 


■ 0.037 


{2a) 


this paper 


Fp(nHc),t2 / 


0.000 ^ 


0.011 - 


■ 0.020 


{2a) 


this paper 


yp(+Hc),t2 ^ 


0.000 ^ 


0.012 - 


■ 0.027 


{2a) 


this paper 



^Recommended values, see text. 
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Fig. 1. — Te(He ii)/Te(0 iii) as a function of Te(0 iii) and temperature fluctuations for the case in 
which all the O is O^^. When O^ is present higher t^ values are expected, particularly for those 
objects with the highest Te(0 iii) values (see Figure 2). Typical t^ values in H ii regions are in the 
0.01 to 0.04 range. 

Fig. 2. — Te(0 III) versus Te(He li) showing the effect of different O"*" fractions and different total t^ 
on Te(He ii). Notice that when O"*" is present there are forbidden regions that correspond to the low 
end of t^ (those regions where the lines stop) because it is unphysical to assume a uniform Te(He ii) 
when one is already using 2 different temperatures for a given line of sight: Te(0 iii) and Te(0 ii). 
The Te(0 ii) for the O^ fraction is assumed to be the one derived from the photoionization models 
by Stasihska (1990). 

Fig. 3. — y(nHc) versus 0/H diagram. Results for t^ = 0.000 (Izotov &: Thuan 1998; Izotov 
et al. 1999) are shown as open circles, the error bars are those quoted by the authors, they include 
the errors in the measurements of the lines, but assume no uncertainties on the determination of 
Te(He ii), A''e(He ii), or r(3889); results for t^ ^ 0.000 (this paper) are shown as solid dots, the 
errors represent the uncertainties on: the intensities of the lines, re(He ii), A^e(He il), and r(3889); 
other possible systematic effects (such as the coUisional contribution to the Balmer lines, or errors 
on the atomic physics) are not included, see text. The lines represent the best fits to the data 
assuming a slope of AY/AO = 3.5 ± 0.9, see text. 

Fig. 4. — Helium, deuterium, and lithium abundances predicted by standard Big Bang nucleosyn- 
thesis computations with three light neutrino species for different values of ry, the baryon to photon 
ratio (Thomas et al. 1994; Fiorentini et al. 1998). Also in this figure we present observational 
abundance determinations of these elements. The helium boxes are la determinations from this 
work; the solid line boxes correspond to l^(nHc) and the dotted line boxes to 1^(+Hc) (see Tables 
9 and 10). The deuterium box is the la determination by O'Meara et al. (2001) while the lithium 
box is the 2a determination by Suzuki, Yoshii, & Beers (2000). 
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